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The existence of conduction channels in lithium sihcate {Li20){Si02) is investigated. Regions 
of the system where many different ions pass by form channels and are thus spatially correlated. 
For a closer analysis the properties of the individual ionic sites are elucidated. The mobility of ions 
in single sites is found to depend strongly on the number of bridging oxygens in the coordination 
shell. The channels are not reflected in the network structure as obtained from the distribution of 
the bridging oxygens. Spatial correlations similar to those found in the silicate also emerge from 
studying the dynamics of particles in a simple random lattice model. This supports the suggestion 
that the observed spatial correlations can be viewed in analogy to the emergence of percolation 
paths. 



I. INTRODUCTION 

The presence of conduction pathways has been pro- 
posed to rationalize the fast conduction in ion conductors 
hke alkah sihcates. Based on EXAFS data on the coordi- 
nation environments in these systems, Greaves proposed 
the modified random network model for glasses^. It pre- 
dicts the aggregation of modifier ions into channels lined 
by non-bridging oxygen atoms. Later NMR studies2i^ 
reported both inhomogeneous and homogeneous distri- 
butions of cations in the structure with a clustering of 
alkali ions in the low concentration regime. These results 
are backed by molecular dynamics (MD) simulations of 
alkali silicate where also only for low alkali concentrations 
a mild clustering of alkali ions has been found^'^. 

In his analysis of MD simulations of sodium silicate, 
Jund et al^ focussed on the most mobile cations. They 
introduced ndiff, the number of different ions visiting a 
given subvolume, to determine those parts of the system 
taking part in effective long range transport. They report 
that the regions of the silicate network, visited by the 
largest number of different ions form a network of blobs 
and connecting channels (denoted conduction channels). 
Stated differently, these regions show a significant spatial 
correlation beyond the nearest-neighbor shell. In particu- 
lar, these spatial correlations are much stronger than any 
possible correlation of alkali positions themselves. Thus 
these conduction channels do not simply emerge from 
clustering of alkalis^ . 

Interestingly, all static structure factors for sodium- 
silicate exhibit a prepeak at qi w 0.95 A^^, correspond- 
ing to the next-nearest Na-Na or Si-Na neighbors, as 
shown experimentalljii and numericallj&S. Thus the 
structure displays some long-range correlations. One 
may be tempted to relate these long-range correlations to 
the channels discussed above. The validity of this conclu- 
sion, however, is not obvious, since the structure refers to 
all ions whereas the channels emerge from the properties 
of the specific subset of regions, visited by many different 
ions. 

In this paper we analyze the question to which de- 
gree these conduction channels are predetermined by the 
structure. This question is central because it elucidates 



the role of the network for long-range transport. Two 
extreme scenarios are conceivable. First, following the 
concepts of Greaves, the network may supply channels 
where ion conduction is strongly favored due to the ag- 
gregation of the relevant structural elements (in Greaves 
work: non-bridging oxygens) which favor fast dynam- 
ics. In this scenario these relevant structural elements 
would also display long-range spatial correlations. Sec- 
ond, the structure might only determine the local mobil- 
ity of an ion such that the relevant structural elements 
do not display these long-range correlations. Rather the 
formation of conduction channels is a simple statistical 
process in analogy to the formation of percolation paths 
in random-energy or random-barrier lattice models of ion 
conductiortiSi. 

In previous work we have presented a method how to 
identify individual lithium sites which are basically time- 
invariant below the glass transitionii. In particular, we 
could show that the continuous trajectories of the lithium 
ions can be interpreted as hops between these sites. With 
this information at hand we can identify the conduction 
channels by counting the number of different ions visit- 
ing a specific site during a long MD simulation. This 
discretization of the system is somewhat more adapted 
to the conduction process than the simple tiling by small 
cubes as done by Jund et al&. In particular, we are very 
sensitive to the effect of structural elements like bridging 
or non-bridging oxygens on the alkali dynamics. 

The paper is organized as follows. In Section ^] we 
discuss the technical aspects of our work. Section IIIII 
contains the results which are finally discussed in Section 

nvi 



II. TECHNICAL ASPECTS 

The data for this investigation are based on MD sim- 
ulations of lithium silicate {Li20){Si02) ■ The system 
containing 1152 atoms was propagated under NVT con- 
ditions at 640K with 2 fs stepsize for 10 ns simulation 
time. As a basis for investigations of the cation dynam- 
ics we first located the individual sites available to the 
lithium ions. For this purpose we have discretized the 
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system in very small cubes (0.3 A) ^ and determined the 
occupation number in every cube. This is a simple way to 
sample the effective potential energy landscape provided 
by the silicate network. In order to distinguish regions 
which indeed belong to ionic sites and those which only 
serve as transition paths between these sites we have used 
some reasonable cutoff criterion for the identification of 
those cubes which are part of ionic sites. Subsequent 
cluster analysis of the remaining cubes finally defines the 
ionic sites. For our system we have obained 378 different 
sites. Details about the method and about the properties 
of the sites are given in an earlier paper—. 

We characterize the mobility of a site in two differ- 
ent ways. First, in analogy to the work by Jund et al. 
the number of different ions visiting a site has been de- 
termined. It is denoted ndig. Sites for which the value 
of ridiff belongs to the 10% highest are denoted A-sites. 
In contrast, those with the 10% smallest values of ndiff 
are denoted A-sites. Second, we determined the average 
residence time in a site, r. Analoguosly, sites with the 
10% shortest residence times are denoted B-sites, those 
with the longest 10% residence times _B-sites. In general 
terms, A- and B-sites indicate positions in the network 
where ions are fast. A- and i?-sites where ions are slow. 
In this investigation both quantities will be used together 
as different measures of site mobility. 

Furthermore we performed straightforward Monte 
Carlo simulations for hopping particles on a 2D square 
lattice model with 1600 sites and periodic boundary con- 
ditions. We introduced random site energy disorder and 
additional barrier disorder among adjacent sites. Both 
distributions were chosen as constant with energies be- 
tween and 8 ksT. 90% of all sites were populated and 
the hopping particles interacted via simple excluded vol- 
ume interaction. The energy disorder was kept constant 
during the whole run. 
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FIG. 1: Scatter plot of the number of different ions ridis, 
vising a site vs. the average residence time of a site r. 
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FIG. 2: Mixed g(r) for Li and network species 



III. RESULTS 
A. General properties of sites 

The two definitions of site mobility are compared in 
Fig. n which shows rtdiff' plotted against r for all sites. 
The highest possible value of ridiff for a site is basically 
given by the number of residences permitted by their 
average duration t: ridis <= t/r. This implies that a 
sufficiently short average residence time is a necessary 
condition for a large number of different ions visiting a 
site. A broad range of lower values for Tidiff is also found, 
caused by multiple visits of identical ions. The correla- 
tion between r and n^iff is therefore only weak. 

General information on the structural environment of 
the sites can be obtained from the radial distribution 
function g(r) of network atoms around lithium. It is 
shown in Fig. |21for bridging oxygens (BOs), non-bridging 
oxygens (NBOs), and silicon atoms. BOs and NBOs have 
been distinguished according to the number of silicon 



atoms within 2.4 A, the distance of the first minimum 
in the g(r) of silicon and oxygen. One finds that lithium 
atoms are surrounded by oxygens. Integration up to the 
first minimum gives average coordination numbers of 4.0 
NBOs and 1.2 BOs. In contrast, if the average number of 
5.2 oxygen neighbors were divided according to the rela- 
tive numbers of BOs and NBOs in the system one would 
obtain 3.3 NBOs and 1.8 BOs. As expected, NBOs are 
favored as coordination partners of lithiumiiSiiii^. 



B. Spatial correlations of sites 

Now we analyze the question to which degree the mo- 
bilities of the different sites are correlated. The presence 
of channels, defined as regions of the system where many 
different ions pass would imply that A-sites (and corre- 
spondingly A-sites) are spatially correlated over larger 
distances. Furthermore, also possible spatial correlations 
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FIG. 3: Correlation of mobilities to the next neighbor shell of 
sites 
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FIG. 4: Correlation of mobilities to the second next neighbor 
shell of sites 



of B-sites are analyzed. 

It can be expected that correlations on the length scale 
of the nearest-neighbor distance occur due to trivial rea- 
sons. For a fast B-site it can be expected that there 
is a low barrier to an adjacent site. Accordingly, also 
the residence time in this adjacent site can be expected 
to be small because a fast backjump to the original site 
would be possible. This trivial correlation is limited to 
next neighbors. If B-sites were arranged into channels, 
they must show an additional correlation over longer dis- 
tances. The presence of a continuous pathway of B-sites 
through the system would result in an increased proba- 
bility to find other B-sites among the second neighbors 
of a given B-site. 

Apart from A- and A-sites (and analogously B- and 
_B-sites) we also selected a subset of 10% of sites with 
intermediate r and ridiff, respectively. Each of the three 
groups contained 37 sites, i.e. ca. 10% of the total num- 
ber. Because of the low number of sites we did not com- 
pute their complete pair distribution functions g(r). In- 
stead we counted directly the number of neighbors from 
each of the three groups around a given central site. The 
limits of the first and second neighbor shell were taken 
from the minima in the g(r) of all sites at 4.2 A and 6.8 A. 
The results for slow, intermediate and fast central sites 
are compared in Fig. O and Fig. ^ for the first and for 
the second shell respectively. 

All values are given as a fraction of the total number 
of neighbors. The statistical value of 37/378 = 9.8% for 
the contribution of each group is shown by the dashed 
lines. Generally, there is a clear tendency that a site 
will be surrounded by first neighbors of similar mobil- 
ity. The strongest correlation is found among the fast 
sites. While the findings are similar for both r and ridifi, 
they are more pronounced in the latter case. For the 
second shell of neighbors, the correlations almost vanish 



for the residence time t, i.e. B-sites (and B-sites) are 
spatially not correlated. In contrast, as already reported 
by Jund et al. A-sites are significantly correlated beyond 
the nearest-neighbor shell. The correlation for A-sites is 
much weaker, but still present. 



C. Spatial correlations of relevant network 
properties 

Next we turn to the question in which sense the con- 
duction channels are predetermined by the structure. To 
answer this question we have to proceed in two steps. 
First, we have to check which property of the network 
significantly favors highly mobile sites. As will be shown 
below it is the presence of bridging oxygens (BOs) which 
displays the strongest correlation with the mobility of 
a nearby site. Second, possible long-range spatial cor- 
relations of this network property (here: occurrence of 
BOs) have to be analyzed. More specifically we check 
whether sites with a high/low number of BOs are spa- 
tially correlated. Again we have to be aware that trivial 
nearest- neighbor correlations will be present. If some site 
displays an increased number of BOs as compared to the 
average site it can be expected that a small number of 
neighbor sites shows the same correlation because one 
BO typically belongs to 2 lithium ions. 

Using the sets of sites already defined above, we first 
count the number of BOs, NBOs and Si atoms up to the 
1st minimum distances of 2.9A, 3.0A, and 4.1 A respec- 
tively. The results are shown in Fig. [S] As expected 
from the above-mentioned analysis of coordination num- 
bers, sites are on average surrounded by « 4 NBOs, while 
only 0.5 to 1.5 BOs are found. For comparison, lines are 
shown for 1.8 BOs and 3.3 NBOs, i.e. the statistical val- 
ues (see above). Interestingly, the strongest infiuence on 
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FIG. 5: Network species as next neighbors of sites 



the mobility of a site is shown for the number of coordi- 
nated BOs, which decreases by « 1 from slow to fast sites. 
The number of silicon neighbors shows only a weak de- 
crease in the same direction. For NBOs the correlation is 
reverse, but much weaker even on an absolute scale. The 
results are very similar when using t or ndiff to determine 
the mobility of the sites. 

This result implies that in most cases a single BO in 
the nearest-neighbor shell is sufficient to immobilize ions 
on the sites around the BO. A possible underlying struc- 
ture in the network that could determine the layout of 
conduction channels should therefore be related to the 
distribution of BOs. Similarly to the method described 
above, we sorted all sites into three groups according to 
the number of BOs in their 1st shell closer than 2.9 A. 
Then the fraction of neighbors from each of these groups 
was counted. The results are given in Fig. The dashed 
line again marks the statistical value, which is 0.33 in 
this case, corresponding to groups of one third of all sites. 
The expected trivial correlation among nearest neighbors 
is, of course, found. Most important, no correlation to 
the second nearest neighbor shell is found. We checked 
that also the NBOs do not display any correlations be- 
yond the nearest neighbor shell. The lack of spatial cor- 
relations beyond the nearest neighbors strongly suggests 
that even the distribution of bridging oxygen atoms as 
the most relevant network property does not determine 
complete conduction channels in the system. 



IV. DISCUSSION 

The presence of only one BO drastically slows down the 
ionic transport through nearby sites. An indication how 
this slowdown is effected can be obtained from the two 
examples shown in Fig[7| Here the local network struc- 
ture is shown around one slow site (left) and one fast site. 
Cations are shown at the center as solid balls. The sur- 
rounding silicate network is shown as tubes, . The site 
as determined by ua^i is represented by the translucent 



0.5 



Oh 
O 

^0.4 

01 

O 



Oh 



0.2 





1 


1 


▲ many 








■ mid 


- 






▼ few (1st sh) ■ 








/\ many 




▼ 




□ mid 








Y few (2nd sh) 








▲ 






A 

A 




--- 


V 


--5-- 


^— -- 








▼ 




▲ 

1 


1 


1 




few 


mid 


many 



group at center 



FIG. 6: Correlations in the number of neighboring BOs of 
sites 




FIG. 7: Local structure around a slow site (left) and a fast 
site (created with the software vmd— ) 



cloud surrounding the cation. Connections are drawn be- 
tween the lithium ion and all oxygen atoms within the 
first neighbor shell. For the fast site, all six coordinating 
oxygen atoms are NBOs, which are bound to four dif- 
ferent silicon atoms. For the slow site there is one BO 
among the neighbors. This BO is part of a silica chain 
segment that is also carrying three of the NBOs coordi- 
nating the cation. These NBOs are bound to three silicon 
atoms on different sides of the BO, with an additional BO 
inserted on one side in between. In total, a chain segment 
of 5 atoms is linked to the cation via four coordinating 
oxygen atoms. The chain is thus closely wrapped around 
a large part of the site. This situation suggests that BOs 
introduce high transport barriers for nearby sites because 
the silica chains joined by them block a large fraction of 
possible jump pathways. 

The pathways through the system used for effective 
ionic transport follow continuous channels of sites visited 
by many ions. The passage of a large number of different 
ions through a specific site depends on two conditions, (i) 
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FIG. 8: Spatial correlations in a simple lattice model 

A comparatively low average residence time of the sites is 
necessary but not sufficient to enable the passage of many 
ions, (ii) The presence of coordinating BOs strongly de- 
creases the number of ions reaching a site, probably by 
introducing a local barrier. Both the distribution of sites 
with low r and that of BOs in the structure show no 
long range order comparable to that found in the path- 
ways traced by sites with high ndig. They can therefore 
only impose constraints on the layout of these pathways, 
but not determine them. Therefore the emergence of 
conduction channels does not seem to require preformed 
channels provided by the network structure. 

Now we show that on a qualitative level similar effects 
can be observed for the simple lattice model, introduced 
in Section ^ After performing long Monte Carlo sim- 
ulations we determined the A- and B-sites similarly to 
our procedure for the silicate system. They are marked 
by filled dots and open circles, respectively. Interest- 



ingly, the three main observations from the silicate sys- 
tem are reproduced: (i) By definition the structure (here: 
the local energetic parameters) does not display correla- 
tions beyond the nearest neighbor shell, (ii) The B-sites 
are not spatially correlated beyond the nearest-neighbor 
shell. Actually, one clearly sees that sites with short 
residence time typically occur as pairs, because of the 
trivial correlations discussed above, (iii) The A-sites are 
spatially correlated, i.e. they are concentrated in a few 
larger blobs. These regions contain some of the groups of 
sites with low r, but also some sites with higher residence 
times. Correspondingly, some sites with low r are not in- 
corporated into regions with high ridiff. Sites visited by 
many different ions can of course not be completely iso- 
lated. The ions have to move away through other sites 
in the neighborhood to allow the following ions to enter. 
In our lattice model the regions best suited for this pro- 
cess are located by chance. The ions benefit from sites 
with very low residence times, they must reach them via 
other sites that are also not too slow, and they must 
avoid high barriers distributed through the lattice. In- 
terestingly, the pathways are somewhat extended in both 
dimensions. This suggests that regions where many dif- 
ferent particles may diffuse are extended in all directions. 
Actually, the channels, as obtained from the molecular 
dynamics simulations^ show this property as well. 

In summary, we have analyzed the nature of the con- 
duction pathways in lithium silicate as well as a simple 
random lattice model. In both cases we do find spatial 
correlations in regions of high mobility. These structures 
do not seem to be supplemented by corresponding spa- 
tial correlations of network properties, following from the 
lack of long-range correlations of sites with BOs. Thus 
we would like to suggest that the formation of conduc- 
tion channels is rather a statistical selection process in 
analogy to the formation of percolation paths in simple 
lattice models. 
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